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Abstract 



We show that the stacking of flat aromatic molecules on a curved surface results in topological 
defects. We consider, as an example, spherical vesicles, self-assembled from molecules with 5- 
and 6-thiophcnc cores. We predict that the symmetry of the molecules influences the number 
' of topological defects and the resulting equilibrium shape. 



Engineering of molecular assemblies is one of the aims of modern nanotechnology. Here, starting 
from ideas put forward by Nelson, MacKintosh and Lubensky for liquid crystals [U [2] [3] we show 
t— ( , that topological defects can determine the shape of supramolecular vesicles [1]. Recently, Nelson [I] 
proposed an elegant way to generate topological defects by coating spherical colloidal particle with 
anisotropic object like nematic liquid crystals. The high-energy core of defects can be used as sites 
for chemical activity and as potential spots for the assembly of three dimensional architectures. 
Further, the classical works [3] relate the appearance of topological defects in lipid bilayers to 
the temperature-induced phase transition from the Sm-A phase to the Sm-C phase with tilted 
^ ■ molecules. 

In the last decades supramolecular chemistry has created large aromatic molecules that form a a 
variety of nanometer sized structures (5] [6j. Here we demonstrate that packing of flat aromatic 
q \ molecules on a curved surface, e.g. on a sphere, results in unavoidable topological defects, similar 
to those considered by Nelson. Unlike colloidal particles and two dimensional spherical crystals 
j— i , [8], the defects in aromatic self-assemblies could influence the geometry and the topology of the 
structure as a whole. Therefore, we allow spherical vesicles to deform towards spheroids in order 
' to determine the minimum energy configuration. 

As an example, consider the bolaamphiphilic sexithiophene (6T) molecules [6j |9j known to self- 
assemble in propanol-2 into spherical vesicles with radius of the order of 100 nm. The aromatic 6T 
core can be thought of as a flat plaquette^ (see Fig. [T^,b). Consider now the regular arrangement 
of flat plaquettes on a curved surface, e.g. a spherical cap (see Fig. Hfc): every plaquette can be 
■ described by a director (headless vector) normal to the plane of the plaquette and thus tangent 
J> , to a sphere. The hairy ball theorem (or generalized Poincare-Hopf theorem) states that there is 
°k> \ no nonvanishing continuous tangent (vector) field on a sphere. This means that we cannot pack 
flat molecules on a sphere without creating a topological defect, namely, a point around which 



the tangent field rotates by 2-7rg, with q an integer or half-integer number, like the one shown in 
Fig. QJi. The number of topological defects depends on the Euler characteristic of the surface and 
on the intrinsic symmetry of the self-assembled molecules. The 6T core, enclosed in a rectangle, 
has 2/m lattice symmetry [ID], which means that there are no invariant vectors in the plane of the 
molecule or that there is inversion symmetry. This defines a unique normal n to the plane of the 
plaquette which will be a tangent vector to a sphere. Instead, five-thiophene groups (5T) have 
mm2 point group transformations, which leave a vector v = (x, 0, 0) invariant under symmetry 
operations. This additional symmetry of the core yields a director to the plane of the molecule 
with n — > — n. We may generalize our results to even and odd number of thiophene and benzene 
rings. For any aromatic molecule, it is always possible to assign either a true vector normal to the 
plaquette or a director. For the tangent vector field on a sphere, there are two topological defects 
with vorticity 27r (topological charge q = +1), whereas for directors there are four defects with 
vorticity ir (q = +1/2) [U [2] - Based on the example of 6T molecules, which can be generalized to 



1 Bolaamphiphile 2,5"" '-(R-2-methyl-3,6,9,12,15-pentaoxahexadecyl ester) 6T has two polar ethylene oxide chains, 
which do not play any role in our consideration. 




Figure 1: (a) Aromatic molecules with six thiophene (6T) cores and five thiophene cores (5T) that 
can be thought of as flat plaquettes (b) Since for 6T there is no invariant vector v, the normal 
to the plane is a true vector n in the z-direction. For 5T we can identify the invariant vector v 
(in x-direction), and thus there is a director n —* — n in the z-direction. (c) The arrangement of 
flat plaquettes on a spherical cap, resulting in a topological defect at the center, (d) The tangent 
vector field and two topological defects (thick dots) located at the north (N) and the south (S) 
poles of a spheroid; c and a are the semiaxes of the spheroid. 



any aromatic molecule with inversion symmetry, we consider in the following the interaction of two 
topological defects on surfaces with the topology of a sphere. 

Like charged particles, the defects interact with each other and with the Gaussian curvature as a 
background charge distribution [111 [T2] . Taking into account the fact that self-assembled structures 
can adjust their shape in order to minimize the free energy, we allow deformations of the sphere 
towards spheroids. This feature is essentially different from the case of colloidal particles where the 
geometry is determined by the substrate pQ . The quantitative understanding of the energetics of 
topological defects interacting on spheroids involves the calculation of the Green function r(xj,Xj), 
which is the inverse of the Laplace-Beltrami operator on spheroids [U [?]. This is a non-trivial 
problem, which can be solved through the conformal mapping of the complex surface onto a simple 
one, e.g. plane or sphere, where the Green function is known |12j . 

To quantify the effect of topological defects we calculate the total free energy F to tai as the sum of 
bending energy Fb en d and of the energy associated with topological defects infect on spheroids, 
namely 

-Ftotal = ^bend + -^defect = 2fc JJdSH 2 + 4lT 2 K A qiqjT (xi,Xj) , (1) 

where H is the mean curvature, k is the bending rigidity, and Ka is hexatic constant [2j 113] . 
For every point of the surface one can define the two principal curvatures k% and «2- The mean 
curvature H = (k% + k%)/2 at a given point depends on the change of the normal vector around 
this point, characterizing the local properties of the surface. The Gaussian curvature K = K1K2, 
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Figure 2: (a) The normalized free energy i^otai/ (^k) as function of c/a (the axes of the spheroid, see 
Fig. [TH) in logarithmic x-scale. The surface area of the vesicle is kept constant, assuming constant 
density of the molecules. The minimum of i^otab corresponding to a sphere c = a for Ka = 0, shifts 
towards higher values of c/a with hexatic constant Ka compared to k. For the ratio KA/k = 1.2 
we expect an instability towards rod-like structures, (b) The equilibrium shapes marked with the 
thick dots on the plot. 



instead, is an intrinsic property of the surface defined only by its metric, and its integral over a 
smooth surface J dS K is related by Gauss-Bonnet theorem [14] to the Euler characteristic x, which 
determines the topology of the surface. Since we consider spheroids with x = 2, the integral over 
K is constant and therefore it is not included in Eq. (pQ). The non-local (long-ranged) interactions 
of the defects are described by the Green function r(xj,x J ), which depends only on the metric 
of the surface. The competition between two terms in Eq. ([I]) leads to frustration, because the 
defects favour high Gaussian curvature [TT] , which inevitably leads to high mean curvature (due to 
H 2 ^ K), and consequently to an increase of the bending energy Fb en d- Because of the repulsive 
nature of defect interactions, the minimum of the second term corresponds to defects located at the 
north (N) and the south (S) poles of a spheroid for a tangent vector field (see Fig. [TH), and at the 
vertices of a tetrahedron for a director field [2]. The topological defects in self-assemblies account 
for the competition between local interactions of aromatic cores and the topological constraints, 
whereas the details of microscopic potentials reside in the bending rigidity k and the hexatic 
constant Ka- In our system, the ratio KA/k is the only governing parameter, related to details 
of intermolecular interactions for a given temperature. Although the value of the bending rigidity 
k = (2.56 ± 0.8) • 10~ 21 J was found from experimental measurements of magnetic deformations of 
6T spherical vesicles in isopropanol [15], the value of the hexatic constant Ka is a less accessible 
parameter. Assuming that tt-tt interactions between thiophene molecules give the main contribution 
to the values of k and Ka we can roughly estimate the ratio KA/k based on the results of ah initio 
molecular orbitals calculations for different geometries of thiophene dimer |16[ I17j. yielding KA/k 
of the order of unity. 

It turns out, that for the case of spheroid the expression for the total free energy F tota j (Eq. (TjQ)) 
can be given in simple analytical form as the function of c/a, the ratio of the semiaxes of spheroids, 



with bending energy 



^bend = 2fc7r(c/a) 2 



tanh -1 y/l - (c/a) 2 



+ 



2 1 



^( 1+ (^f)}' (2) 



(c/a) 2 3(c/a) 



and the Green function 



r(o,7r) 




log — + \J -\ + (c/a) 2 tan 



V-l + (c/a) 2 



(3) 



Note that since the Gaussian curvature at the poles is given by K = c 2 /a 4 , T(0, 7r) becomes more 
negative for increasing K. We plot in Fig. [2^ the total free energy i^totai for different values of 
Kji/k as a function of c/a with defects located at z = c (N) and z = — c (S). In presence of 
topological defects, the equilibrium shapes, corresponding to the minimum of -Ftotah are shown for 
two values of KA/k (see Fig. [2b). When the contribution infect is not dominant (Ka % 0.6/s) 
the equilibrium shapes are close to spherical, but for K a — k they are elongated with c/a ~ 3.2. 
For small values of the bending rigidity (Ka > k) we expect an instability towards infinitely long 
tubular structures. A similar elongation of phospholipid vesicles has been predicted to occur as a 
function of temperature at the phase transition between Sm-A phase and Sm-C phase [3]. Indeed, 
6T molecules in butanol, in contrast with isopropanol, were found to self-assemble into multiwalled 
cylindrycal structures [H]. The theory presented here holds only for single layered structures and 
thus we cannot establish a direct comparison with the experimental observations. Moreover, in real 
situations, parameters such as k depend also on the particular solvent |19j . 

In conclusion, we demonstrated that the stacking of flat aromatic molecules in self-assembled spher- 
ical vesicles leads to topological defects. The essential difference in symmetry between molecules 
with even and odd number of thiophene groups results in two and four topological defects, re- 
spectively. In the case of two interacting defects the equilibrium shape of vesicles turns out to 
be elongated. We believe that four defects at the vertices of a tetrahedron would hardly lead to 
an elongation, although a detailed analysis of four interacting defects is beyond the scope of this 
paper. In general, the question of how the shape of self-assembled structures is connected with 
the intrinsic properties, such as the symmetry, of the constituent molecules remains open. Here we 
proposed a simple geometric approach based on topological defects to establish this connection for 
molecules with thiophene/benzene aromatic cores. The proposed theory can be used to extract the 
value of KA/k from experimental observations of the shape of self-assembled aromatic molecules. 
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1 Point groups for 6T and 5T 

In Figure la we show two aromatic molecules (6T and 5T) with even and odd number of thiophene 
cores and the corresponding symmetry operations, which transform the molecule into itself. For 
two-dimensional rectangular lattice, like the one enclosing 6T and 5T molecules, the possible sym- 
metry operations are mirror-reflections across line, denoted by to. For 6T molecule, we can write 
the following symmetry operations which leave one point unmoved [10] 
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m z =| 1 j and 2 2 = m x m y = ^ —1 | . (4) 

Their product m x m y m z = I —1 | , acting on a vector, results in the total inversion 
symmetry. For 5T molecule the point group is different, containing the operators 

A 

, 2 X = m y m z = I -1 | (5) 

These symmetry operations leave vector v = (x, 0, 0) invariant, which is not the case for 6T 
molecules. This result can be generalized to even and odd number of thiophene and benzene 
rings. 

2 Conformal mapping 

Let us consider a spheroid, which is a surface of revolution given by the following parametrization 

x(u,v) = (a sin v cos u, a sin?; sin??, ccosv), < ?; < 7r, <u < 2ir, (6) 

with semiaxes c and a (see Figure Id). Then, the metric of spheroid is given by 

ds\ = (a 2 cos 2 v + (? sin 2 ?;) dv 2 + a 2 sin 2 ?; du 2 . (7) 

The mapping of spheroids onto a sphere with metric ds 2 = R 2 (d9 2 + s'm 2 6d(j) 2 ) is said to be 
conformal if we can write 

ds 2 = e 2X ^ds 2 , (8) 

where e 2A ( u ** is called the conformal factor, which varies with position u = {u, v} on spheroid. 
Because of the rotational symmetry u = (j) an d conformal factor depends only on parameter v. By 
equating two metrics we find that 



o\( v \ R 2 sin 2 d6 , , / 9 c 2 , n . 

e 2X{v > = , =±dv\ cot 2 v + ^. 9 

a z sin v sm a V a z 



By integrating both sides of the second equality we get 



log ( tan — ^ = g(v) = — \J — 1 + rj 2 arctan [ = - | — log ( cot v + a/ 7 / 2 + c °t 2 v ) , (10) 

V 2 / \ V rf + cot 2 v J 



rj = c/a 



where 



expression for t 



Substituting sinfl = 2tan(0/2)/(l + tan 2 6/2) into Eq. © we find the analytical 
re conformal factor 

e 2A W = — , (11) 

a 2 sivPv cosh 2 <7(t>) ' 



and the limit of interest 



R 



lim X(v) = log \J — 1 + rf arctan a/— 1 + rj 2 . (12) 

d^o a 

3 Free energy 

3.1 Contribution from defects 

The pair Green function for topological defects on deformed sphere can be defined similar as for 
superfluids in reference [12] 

r( Xi , Xj ) = -i-bg^ + -^(Afo) + A( Xi )), (13) 

where T>ij is the distance between two defects on the sphere (the chord between two points!) and 
d is the core size of the defedH. This representation of the Green function already includes both 
the interaction between defects (first term) and the position dependent self-energy of defect. In 
the case of two topological defects on sphere with topological charge g, = g,- = +1, the minimum 
of -Fd e f ec t = Air 2 K Aqiqj^{x.i,~x.j) corresponds to the defects located at the north N (v = 0) and the 
south S (v = 7r) poles (see Figure Id). It gives = 2R, and together with Eqs. (fl"2|) and (fT3j) we 
find the Green function 



r(0, 7r) = ~^:^ lo S~^ + + V 2 arctan \J —I + rf^j , 



(14) 



and consequently an analytical expression for infect- In Figure [3] below we plot the Green function 
r(0, 7r) and the Gaussian curvature K as the function of rj. As was expected, the Green function 
decreases with rj, while the Gaussian curvature is increasing, resulting in the more negative r(0, tt) 
and thus -Fdefect f° r higher K, which is in agreement with a theory |llj . 

3.2 Bending energy 

The bending energy proposed by Helfrich in |13| is written as 

F head = 2k [fdSH 2 , (15) 



where H is the mean curvature and the integral is over the surface S. For spheroids (Eq. [6]) 



dS = a sin v\J a 2 cos 2 v + c 2 sin 2 v du dv, (16) 

h=-( ^ + — y (17) 

2 I y/ a 2 cos 2 v + c 2 sin 2 ~ (a 2 cos 2 v + c 2 sin 2 vfl 2 J 



(18) 



After some calculations, the integral in Eq. (115j) can be simplified to the following form 

, 2 f tanh" 1 Jl-7] 2 2 1 / 2\1 

2 At this microscopic scale the continuum theory breaks, so d may be thought as a cut-off radius. 



(a) T(0, it) 



(b) K 



-2 
-3 
-4 

-5; 

-6 



150 
100 
50 



0.2 0.5 1.0 2.0 5.0 10.0 



0.2 0.5 1.0 2.0 5.0 10.0 



V 



Figure 3: (a) The Green function T(0, ir) given by Eq. (|14|). For this plot we choose the equilibrium 
radius of the vesicle R = 100 nm and the size of the core d = 3.5 Aas the equilibrium distance 
between thiophene cores, (b) The Gaussian curvature K = c 2 /a 4 = rj 2 /a 2 calculated at the poles, 
assuming the condition of constant surface (Eq. [T9|) . 



This energy term does not depend on the size of the vesicle only on the dimensionless parameter r\. 
However, the Green function in Eq. (114j) depends on the absolute value of the distance between the 
defects, therefore in all our calculations we assumed the condition of constant surface of spheroid 

5 = 2 7 ra 2 fl + r ? 2tanh " lv/r ^Uconst, (19) 
V V 1 - rf J 

yielding constant density of the molecules. 



4 Estimate of the ratio of KA/k 

We assume that for aromatic molecules this value is mainly determined by the long-ranged ir- 
ir interactions between aromatic rings, and not by extrinsic effect (e.g. solvent) as discussed 
in |17] for thiophene oligomers. Therefore, we propose to estimate the ratio K^jk by using 
the results of quantum chemistry calculations for the binding energy of thiophene dimers with 
different relative orientations [16]. The basic idea is to associate the energy of the splay con- 
figuration of two iV-thiophene molecules with the value of the bending rigidity k (see Fig. [3b 
below), and the value of hexatic constant K~a with the rotation around the long axes of the 
molecule (see Fig. 0b below). Then, based on the calculated interaction energy for different 
geometries of thiophene dimers (configurations A, G, H of Figure 2, Table 2 of reference |16j ) 
we find k oc E totaX (G) - E totaX (A) = -2.05 - (-1.32) = -0.73 kcal/mol m 5 • 10" 21 J and 
K A cc E total (H) - E total {A) = -2.28 - (-1.32) = -0.96 kcal/mol, yielding K A /k w 1.3. 
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